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Abstract 

The dynamics near the top of a potential barrier is studied in the temperature 

region where quantum effects become important. The time evolution of the 

density matrix of a system that deviates initially from equilibrium in the 

vicinity of the barrier top but is in local equilibrium away from the barrier 

top is determined. Explicit results are given for a range of parameters where 

the nonequilibrium state is not affected by anharmonicities of the barrier 

potential except for the barrier height. In particular, for a system confined 

initially to one side of the barrier the relaxation to a quasi-stationary flux state 

is determined. The associated rate constant is evaluated and the relation to 

other rate formulas is discussed in detail. 
PACS numbers:82.20.-w, 05.40.+j 
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I. INTRODUCTION 



The dynamics of systems hindered by a potential barrier plays an important role in almost 
all areas of physics and chemistry. The reaction coordinate which describes the transition 
across the barrier typically interacts with many degrees of freedom. In the classical region, 
i.e. for high temperatures, the generalized Langevin equation for the reaction coordinate 
usually provides an adequate description of the barrier dynamics. Based on these stochastic 
methods, Kramers flux over population approach enables a detailed investigation of the 
escape process across the barrier M. 

A corresponding formulation of escape processes in the presence of quantum mechanical 
effects is available only since recently |H§. In the classical region where the barrier is crossed 
by thermally activated processes only the barrier height and the curvature of the potential 
at the barrier top and the well minimum are relevant for the rate constant In this 
article we study a region where quantum effects lead to large deviations from classical rate 
constants but where the harmonic approximation for the barrier potential is still sufficient to 
determine the dynamics of the nonequilibrium state. This allows for analytical results. We 
extend earlier results on quasi-stationary states of systems with large barriers to include the 
short time dynamics and the relaxation to a quasi-stationary flux state. Furthermore, the 
approach will be used to determine time correlation functions of population or flux operators 
associated with the escape process. This provides the connection with other familiar rate 
formulas. 

The article is organized as follows. In section || we give a brief outline of the formalism 



and collect results which are of relevance in the following. In section |TTJ the time evolution 
of an initial state which is in nonequilibrium near the top of the barrier is investigated. 



In section [TV| these results are applied to determine expectation values, as e.g. the avarage 
flux across the barrier, and the flux-flux correlation function. The results are illustrated by 
an explicit example. Finally, the relation to other approaches of quantum rate theory is 
discussed. 
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II. DYNAMICS NEAR THE BARRIER TOP 



In this section we collect some results on the description of dissipative systems that are 
needed in the following sections and introduce basic notation. 

A. Dynamics of dissipative systems 

The stochastic motion of a classical particle of mass M moving in a potential field V(q) 
coupled to a heat bath environment is described by the generalized Langevin equation 

Mq(t) + M jTdf 7 (* - 09(0 + = £(*)■ (1) 

Here, the stochastic force £(t) and the nonlocal damping kernel j(t — t') are connected by 
the fluctuation-dissipation theorem 

mat')) = kBTM 1 (\t-t'\) (2) 

where T is the temperature of the environment and k-Q denotes the Boltzman constant. In 
this paper we consider systems where V(q) has a smooth potential barrier. Then, near 
the barrier top the barrier potential can be approximated by the potential of an inverted 
harmonic oscillator. Assuming that the barrier top is at q = and V(0) = 0, the barrier 
potential may be written as 

via) = ~Muy. (3) 

Within the range of coordinates where this form of the potential is valid, the classical barrier 
dynamics can be determined exactly by means of the Langevin equation ([J). In particular, 
the dynamics near the barrier top depends on local features of the barrier potential only 
and is not affected by anharmonicities. However, when the temperatures is lowered quantum 
effects become important and the barrier dynamics may depend on global features of the 
potential field. 
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The dynamics of a quantum statistical system is determined by the time evolution of the 
corresponding density matrix. Starting at t — from a general initial state Wq of the entire 
system composed of the Brownian particle and the heat bath, one has 

W{t) = exp(-iHt/h)W exp(tHt/h) (4) 

where H contains the Hamiltonians of the system, the environmental degrees of freedom, 
and a system-environment coupling. We shall assume that the state Wq is out of thermal 
equilibrium due to a preparation affecting the degrees of freedom of the Brownian particle 
only. Since we are interested in the dynamics of the particle only, the time evolution of the 
reduced density matrix pit) = tr R W(t) will be considered, where tr R is the trace over the 
reservoir. To eliminate the environmental degrees of freedom it is convenient to employ the 
path integral approach The environmental degrees of freedom can be integrated out 

exactly if the heat bath consists of harmonic oscillators which are coupled linearely to the 
coordinate of the particle. In the limit of infinitly many bath oscillators with a continous 
frequency spectrum this model causes dissipation and in the classical limit the generalized 
Langevin equation (|I|) is recovered. The details of the path integral representation of the 
reduced density matrix and explicit calculations are given elsewhere [[|. As a result, the 
position representation of the time dependent reduced density matrix is found to read 

P(Qf,q'f,t) = J dqidq-dq dq' J(q f , q f , t, q u q u q , q') A(ft, ft, q , q'). (5) 

Here, J(qf,qf,t,qi,qi,q,q') denotes the propagating function represented as a 3-fold path 
integral where two path integrals are in real time arising from the two time-dependent 
operators in (||) and one in imaginary time describes system-bath correlations in the initial 
state. Since for the parabolic barrier (^) the propagating function is given explicitly below, 
we omit here its general form and refer to ||. Equation ([5]) determines the time evolution 
of the density matrix starting from the initial state 

p(?/>S/>0) = dq dq' \(q f ,q,q' f ,q')pp(q,q'), (6) 
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where pp = tr^Wp). Here, Wp is the equilibrium density matrix of the entire system and 
A(g/, q , q'f, q') is a preparation function describing the deviation from thermal equilibrium. 
In an initial state of the form (||) the system and the bath are correlated. Hence, the 
customary assumption that the initial density matrix Wq factorizes into the density matrix 
of the particle and the canonical density matrix of the unperturbed heat bath is avoided. 
This is a crucial point since (H) allows for the investigation of the dynamics of realistic 
physical systems also for short times where preparation effects are important. 

B. Reduced density matrix for an inverted harmonic oscillator 

In || we have shown that anharmonicities of the barrier potential are always essentiell for 
very low temperatures. Here, we investigate the region of high to intermediate temperatures 
where the parabolic approximation (|3]) for the barrier potential is sufficient but quantum 
effects may be important. 

For the harmonic potential (|3]) the path integrals involved in the propagating function 
can be solved exactly. The explicit calculation is performed in [0]. One finds 



p(x f ,r f ,t) = j dxidridxdr J(xf,r f ,t,Xi,ri,x,r) X(x i: r i: x, r) (7) 
where we have introduced sum and difference coordinates 

x = q — q', r = (q + q')/2 (8) 
for qf, q'f and qi, q[ as well as for q,q', respectively. For the propagating function one obtains 



- -x 1 1 1 M ( ~ 2 

J(x f ,rf,t,Xi,n,x,f) = - — == J— j- H u„u r 



Z4*\A(t)\ ^Kp\A\U^ 2 P Vn=i 



x exp ( ^Ep(x, r) + j^t(xf, r f , t, x i} r h x, r) ) . 



(9) 



Here, 



is the well-known minimal imaginary-time action of a damped inverted harmonic oscillator 
at inverse temperature f3 = l/k&T where 

oo 



a = 4 E u « (ii) 



-oo 



and 



1 oo 

n =f-pT, (WnmWnl) - U 2 ) U n . (12) 



Furthermore, 



27m 



hp 

are Matsubara frequencies and 



(13) 



Un= + K|7(K|) -^o) ■ ( 14 ) 

'y(z) denotes the Laplace transform of the macroscopic damping kernel 7(s) which is deter- 
mined by the spectral density I{uj) of the heat bath 

. 2 f°° dul(u) . . . . 

7(s) = m/„ V1T cosM - (l5) 

We note that for a harmonic oscillator the functions A and Q correspond to the variance of 
the position and of the momentum, respectively. However, for a barrier there is no obvious 
physical meaning since e.g. for high temperatures one has A < 0. When the temperature 
is lowered |A| becomes smaller and vanishes for the first time at a critical temperature T c . 
As seen from (|J) and (|T0|) this leads to a divergence of the propagating function. Hence, as 
already discussed in 0, the harmonic approximation is limited to temperatures above the 
critical temperature T c . For temperatures near and below T c anharmonicities of the barrier 
potential field are always essential 0. 

Apart from the pre-exponential factor the time dependence of the propagating function 
is contained in the second part of the exponent of (§). One finds || 



E t (xf,r f ,t,Xi,ri,x,r) = 



-rXiM 



S 



M 2 



A(t) _ 

' A(t) 2AA(t) j +rXf % 



( ... , i(t) 2 \ 5 sim 



+ixxiM -Q 
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IXXf- 



M 2 



+ -x 2 M 



+iXiXf- 



n 

M 2 



2A(t) J "~~ J h 
S h 2 A 



S(t) 



A(t) 



1 - 



A(t) 
M 2 S(t) 2 \ 



S(t) 



h 



A(t) ' AM 2 A(t) 2 y h 2 A 2 
h 2 A 



S(t) - 4rv^(t) 



+-x 2 M 



Q + A 



A(t) 
A{tf 



2M 2 A{tf 



h 2 A 2 



A(t) 



S(t)S{t)M 2 



M2 ( s(t)-^s(t) 



(16) 



'A(t) 2 h 2 A\""" J A(ty 

Hence, the dynamics at a parabolic barrier is essentially determined by the functions A(t) 
and S(t). They are given by the Laplace transforms of || 

n 



A(z) 



2M 



z 2 + z A f{z) - ujI 



(17) 



and 



2 00 z 
S{Z) = ^z^v 2 



A( Z ) - A(\u n \)) . 



(18) 



Within the harmonic approximation the above formulas (f7|) — ( 18|) determine the time 
evolution of the density matrix near the top of a potential barrier starting from an ini- 
tial state with a deviation from thermal equilibrium described by the preparation function 



X(xi,ri,x,r) 



III. DYNAMICS OF THE ESCAPE PROCESS 

Now, we consider a system in a metastable state which may decay by crossing a potential 
barrier. We imagine that the system starts out from a potential well to the left of the barrier. 
Metastability means that the barrier height V& is much larger than other relevant energy 
scales of the system such as k^T and huo, where TuiJq is the excitation energy in the well of the 
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inverted potential. In the temperature region where anharmonicities can be neglected, i.e. 
for temperatures sufficiently above T c , the time evolution of an initial nonequilibrium state 
near the barrier top can be calculated with the propagating function ([)]). In particular, for a 
system prepared at t = in thermal equilibrium in the metastable well, the relaxation to the 
quasi-stationary state with constant flux across the barrier can be investigated. This will be 
done in this section. The stationary flux state was already determined in || by evaluating 
the propagating function in the large time limit. These investigations are extendend in the 
following to include the short time dynamics and the relaxation to the quasi-stationary 
state. Firstly, in |III A| we introduce the initial preparation. Then, in [JIBl we determine 



the time dependent density matrix, and in [111 C| the relaxation to stationary nonequilibrium 
state is investigated. 

A. Initial preparation 

The initial nonequilibrium state at time t = is described by the preparation function 

B 

X(xi, x, f) = S(xi- x)8(ri - f)e(-ri) (19) 

so that the initial state is a thermal equilibrium state restricted to the left side of the barrier 
only. Then, according to (0), the dynamics is given by 

p(x f ,r f ,t) = J dxidn J(xf,rf,t,Xi,ri) Q(-n) (20) 

with 

J(Xf,Tf,t,Xi)Ti) = J{Xf, Tfj tj X{, Tj, Xj,, 7*j). (21) 



In this case the time dependent part of the exponent in the propagating function (|ST 
simplifies to read 

t t (x f ,r f ,t,Xi,ri) = 5± t {xf,Tf,t,Xi,ri,Xi,Ti) = 

8 



XfTfM + XiTf 



+ -x 2 M 



A(t) 
-Q 
hA 



2A(t) 

h 2 A 



rx MS{t) I r, M ? (f) S{t)A{t) ' 
* l 2AA(t) + 1 f h { A A A{t) 



AM 2 A(t) 2 



%0b %3u J 



+ -x 2 f M 
2 f 



2A(t) 2 

n + A 



A(t) 



'M 2 S{tf 



h 2 A 2 

A(t) 2 M 2 



M 2 S(t) 2 * 
h 2 A 2 

- l) ~ A(t) 



S(t)S(t)M 2 



A(ty 



Ti 2 A 



A 2 h 2 

x 2" 



(22) 



B. Time dependent density matrix 



Since the exponents (|10|) and (p2| ) in the propagating function are bilinear functions of the 
coordinates, the integrals in (p0|) are Gaussian and can be evaluated exactly. For large times 
this calculation is performed in detail in ||. For arbitrary times we may proceed accordingly. 
After determining the extremum of the exponent in the propagating function ( f20|) with 
respect to Xi and r^, one first evaluates the Xj-integral. Then, after simple manipulations of 
the remaining rj-integral, the time dependent density matrix may be written in the form 



p(x f , r f , t) = pp(x f , 77) g(x f , 77, t). 



(23) 



Here, 



pp(x,r) 



1 



M 



(24) 



is the equilibrium density matrix for an inverted harmonic oscillator and 

1 ru{x,r,t) 

g(x,r,t) = —= / dzexpl— z 
yvr J -00 v 

= -erfc [— u(x, r, t)] (25) 
is a form factor describing deviations from equilibrium with 

u(x ' r - t)= fm( l -wm) ' (- r+iW W) x )- (26) 

Clearly, the harmonic approximation is valid only for high enough temperatures. For tem- 
peratures near the critical temperature T c where |A| vanishes, the above result becomes 
divergent. 
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C. Relaxation to stationary nonequilibrium state 



Now, we investigate the dynamics of the density matrix ( p3| ) starting from the initial 
state at t = in greater detail. Note that the time dependence of the form factor (^) is 
completely determined by the function S(t). 

Firstly, let us consider small times uot <C 1. There, one has jSj 

which leads to 

l -iS?=w\ t2+0{t3) - (28) 

Then, the function u(x,r,t), which gives the upper bound of integration in (|25|) , reads 



M 1 Mtt 

Hence, using the asymptotic formula 

[■CO ] 

/ dxexp(— x 2 ) ~ — exp(— z 2 ) for Re{^} — > oo (30) 
J z 2,z 



where Re denotes the real part, the leading order expression for the form factor (|25|) in the 
limit u>ot <C 1 is found to read for finite r 

. ^, . / hn t ( Mr 2 Mxr Mtt 2 \ 
,(*,r,t) = 8(-r) + \/^;exp (-_ + ,_ + (31) 

while for r = 

1 1 Ax-JMQ/lh 

g(x, 0,t) = - + -= dz eM-z 2 ) + 0{t). (32) 

2 y/TT JO 

Clearly, for t — > 0+ and r ^ the form factor reduces to the function contained in the 
initial preparation ( |T9"D as expected. On the other hand, at r = the t — > 0+ limit differs 
from the t — > 0— limit by an imaginary part due to the discontinuity of the G function. 
Defining the width A(t) in position space of the nonequilibrium state ( |2"3"D as that value of 
|g|, q < where w(0, q, t) = 1, one gets 
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This reduces to A(t) = y2hfl/Mt for small times in accordance with ([Sip . 

In we have shown that for large times the time evolution of the density matrix near the 
barrier top has a stationary solution. Here, we regain this result from (|23|). Evaluating the 
functions A(t) and S(t) for times larger than 1/cur one gets to leading order an exponential 
growth H according to 

m = -W, ^HM^fW exp( " R() ' (34) 

and 

S(i) = ~m cot( ^f £) ^r^dr^rW) exp(WRt) ' (35) 

Here, j f (z) denotes the derivative of j(z), and cjr is the Grote-Hynes frequency |§ given 
by the positive solution of u; R + co>r7(co>r) = Uq. Eqs. (0) and (|35D describe the unbounded 
motion at the parabolic barrier with corrections that are exponentially decaying in time (see 
|| for details). Hence, the function u(x,r,t) in (f26|) becomes independent of time 



and the density matrix (^) reduces to the stationary nonequilibrium state derived in 0. 
This time independent state describes a constant flux across the potential barrier and gener- 
alizes the well-known Kramers flux state to the temperature region where quantum effects 
are important. The width A(t) from (^) saturates for large times at the finite value 



Aoo = V^w~ (37) 

which coincides with the width of the diagonal part of the equilibrium distribution (f24|). 

From the above discussion it is obvious that a lower bound of time where the stationary 
flux solution holds derives from u^t ^> 1. For very long times depletion of states inside the 
potential well leads to a flux decreasing in time. Hence, for very long times anharmonicities 
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of the barrier potential become important. For a barrier potential with a quartic term as 
leading order anharmonicity the upper bound of time where the density matrix (E3) is valid 



has been estimated in [[|. One obtains the condition exp(u;Rt) <C q a y2Mu>o/h\A\ where q a 
denotes a characteristic length indicating a typical distance from the barrier top at which 
the anharmonic part of the potential becomes essentiell. 

The density matrix (|2~3D depends on local properties of the metastable potential near 
the barrier top only. On the other hand, the metastable state is assumed to be in thermal 
equilibrium near the well bottom. This means that the solution ( p3l) must reduce to the 



thermal equilibrium state for coordinates qf, q'j on the left side of the barrier at distances 
small compared with q a . Now, for t = the equilibrium state extends to the top of the barrier 
and the matching to the equilibrium state in the well is most critical for the stationary flux 
state where A(t) is largest. However, this latter case was examined in 0. One obtains the 
condition 



a '<<rM«> (38) 



where Vb is the barrier height with respect to the well bottom. From a physical point of view 
(j38[) defines the region where the influence of the heat bath on the escape dynamics is strong 
enough to equilibrate particles on a length scale smaller than the scale where anharmonicities 
becomes important. Only then nonequilibrium effects remain localized in coordinate space to 
the barrier region also for longer times. Especially in the classical region where k^T 3> huj® 
and for Ohmic damping 7(2:) =7 Eq. (|38|) reduces to the well-known Kramers condition |]J 
kBTu /Vb <C 7. Here, 1 — ~ 7 for small damping has been used. When the temperature 
is lowered |A| decreases and the range of damping where the stationary solution ( PB| ) is valid 
becomes larger. This is investigated in detail in [T2J. 

IV. DECAY RATE AND RELATION TO OTHER APPROACHES 



In this section the time dependent density matrix derived above is used to evaluate 
expectation values, in particular the average flux across the barrier. Further, the relation of 
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the theory to other approaches to rate constants is discussed. 



A. Average flux and decay rate 

Clearly, the solution (p3| ) contains all relevant information about the nonequilibrium 
state. Now, we want to evaluate the total probability flux at the barrier top q = 0. One has 

J(t) = ^(p5(q) + 5(q)p)t (39) 

where the expectation value {-) t is calculated with respect to the time dependent nonequi- 
librium state. From (|39| ) one has in coordinate representation 



(40) 

x f =0 



Since the essential contribution to the population in the well comes from the region near 
the well bottom, the normalization constant Z in (|40|) can be approximated by the partition 
function of a damped harmonic oscillator with frequeny uj w at the well bottom, i.e. 

z =A^(fi 2-1 FT? 2 ) «p(W- (41) 

Here, Vj, denotes the barrier height with respect to the well bottom. Note that the potential 
was set to at the barrier top. Inserting (p3| ) for rj = and (|4lD into ( [40"D one obtains 

J(t) = r rj(t) (42) 

where 

T = lim J(t) 

denotes the decay rate of the metastable system in the well. We recall that the Grote- 
Hynes frequency cjr is given by the positive solution of o; R + 0^7 (cur) = Uq. The rate ( f43|) 
describes thermally activated transitions across the barrier where the prefactor takes into 
account quantum corrections B JToyTTH . For the time dependent function r](t) one gets 
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This way we have found an analytical result for the dynamic behavior of the average flux 



which is usually studied numerically, see e.g. |T2[. For long times LU^t 3> 1 the above function 
approaches 1. For very small times one obtains from ( $F7\) 



"W^\/^M +0((! ) < 45) 

which gives a finite flux for t — > 0+ while, according to the initial preparation the limit 
t — ^ — leads to a vanishing flux [see also (|3l| ) and (|32])1. Specifically, for finite damping 



wr V w ol A l 



V(0) = — J-2T7T ( 46 ) 



is always larger than 1. As a consequence, the probability flux for t — > 0+ exceeds the rate 
For very high temperatures where hp <C 1, Eq. ([IB) reduces to 7/(0) = The 
corresponding probability flux J(0) = T/cur coincides with the result of classical transition 
state theory [|TJ 

r cl =^ex P (-/3V 6 ). (47) 

Here, we have used the fact that the term in brackets in the prefactor of (^3|) approaches 1 
for hp 1. For lower temperatures |A| decreases and rj(0) becomes larger than 1/Vr. 



B. Flux flux correlation function 

The propagating function can also be used to determine correlation functions. Here we 
consider the right-left spatial correlation function 

C Rh (t) = tr {e[q{t)]e[-q]p„} = (&[q(t)} Q[-q]) p (48) 

where G(-) denotes the step function. Time derivatives of Crl(^) lead to further correlation 
functions, in particular the flux-flux correlation. Below we will see that these correlations 
are connected with other rate formulas. 
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Now, let us evaluate Crl(£) explicitly. Within the presented real time approach this 
correlation function may formally be looked upon as the expectation value of 9(g) at time 
t of a system with an initial "density matrix" 9(— q)pp. The corresponding preparation 
function then takes the form 

X(xi,ri,x,r) = 6 (-r» - Xi/2) 6(xi - x) 5(r 4 - f). (49) 

This way, using (0), the correlation function may be written as 

Crl(*) = J dr f dxidri 9(r/)6 (-r; - 2^/2) J(0, r/, t, r f ) 

= J dr f dxM 9(r / )9(-r0 J(0, r /; £, x h r[ - Xi /2) (50) 

where the propagating function J(xf,rf,t,Xi,ri) is given in (|2lD . We proceed as in section 



[II B and first evaluate the Xi and afterwards the integration. Here, the maximum of the 



exponent in the propagating function with respect to and r- lies at 

4 = 







M 



h [S(t)+iA(t)]%. (51) 

Introducing shifted coordinates X{ = — x® and f\ = r\ — rf a straightforward calculation 
shows that 

Epfa, r[ - xi/2) + S(0, r f , t, x h r- - Xi/2) = 

W) + - + « - [*(,) + *(«)] . (52, 



8AA(t) 2 ^ WJ M 2 J 2A 2AA(i) 

The Gaussian integrals with respect to and are now readily performed. Finally, after 
some further manipulations, we end up with 

I | / 00 \ roc roc 

Crl(*) = 7^-5 II ^n^n / dxexp(x 2 ) / exp(-y 2 ) 
z 7rfi,/5 y Jo Jx/z(t) 

(53) 



where 
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For t 



one has from fllTD 



A(t) = 



(55) 



Hence, z(t) tends to zero and Cm,(t) vanishes for t —>■ as expected. Now, the time derivative 
of (0) yields 



C RL (t) = (F(t)Q{-q))p 

vt u„ I 

{[S(t) +iA{t)] 2 -h 2 A 2 /M 2 } 



1 1 ( TT 2 \ |g(t)|+»|i(t)| 



where 



F= l -\p5{q) + 5{q)p] (57) 
is the flux operator. Finally, a second time derivative gives the flux-flux correlation 



C RL (t) = (F(t)F) f 



\S(t)\+i\Mt)\ [\S(t)\+i\A(t)\] 2 [S{t)+iA{t)] 
{[S(t) + iA(t)} 2 - ft 2 A 2 /M 2 } V2 {[S(t) + iA{t)Y - h 2 AyM 2 } 3/ ' 

The above three correlations are related to the escape rate out of the metastable well as will 



be seen in section IV D 



C. An example: Drude damping 

To illustrate the above results we now consider a Drude model with 7(i) = 
7Co>d exp(— uj)t) by way of example. Clearly, in the limit ujd ^> u;o,7 the Drude model be- 
haves like an Ohmic model execpt for very short times of order 1/ujy. The Laplace-transform 
of j(t) reads 
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7(z) = 7 



(59) 



Then, from (|TTD and (|12|) we obtain 

oo 



and 



A = J_ ^ 1 

n/3 + |z/ n |( 7 ^ D /u; D + |z/ n |) - ul 



-| oo 



(60) 







|(7^ d /u;d + 






z/ 2 

n 


+ 


^ n (7 w d/wd + 


^n|) - ^0 



(61) 



The time dependence of the nonequilibrium state is completely determined by the func- 
tion S(t) in (HH) . Some of the algebra needed to evaluate S(t) for a Drude model explicitly 



is provided in recent work [13]. We obtain 



S(t) 



M 



E 



| cot f-yM exp(V) 



(62) 



Here, Aj, i = 1, 2, 3 denote the poles of given by the three solutions of 

z 3 + u; D 2 2 + z("fLU D - tu 2 ) - w D = 0. 
For the coefficients Cj one has 

Cl = (XI - A 2 )/0 
c 2 = (XI - A 2 )/0 

v2 



(63) 



c 3 = (Af - A 2 )/0 



(64) 



where 



= (Ai - A 2 )AiA 2 + (A 2 - A 3 )A 2 A 3 + (A 3 - Ai)AiA 3 . 
Further, we have introduced the time dependent function 

7^D 



(65) 



at) 



2 x \u n \ exp(-\u n \t) 



E 



W n=^oo(A 2 -^)(A|-^)(Al-^)- 
which can also be written in terms of hypergeometric functions as 



(66) 
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1 A c. 



F{l±;l + ±;e-*)-F{l,-±;l-±;e-*) 

V V V V 



(67) 



With these results for A, Q, and S(t) and a Drude frequency uj-d = lOOc^o we have in- 
vestigated the time evolution of the nonequilibrium state numerically. In Fig. 1 the width 
A(t) of the nonequilibrium state in position space, given in fl33|), is depicted as a function 
of t for various temperatures. For high temperatures damping effects are relevant for inter- 
mediate times only while for lower temperatures they are essentiell for all times. For small 
times A(t) grows faster for stronger damping and reaches a larger asymptotic value for large 
times. This is due to the quantum mechanical effect that stronger damping suppresses the 
fluctuations of the coordinate and therefore enhances fluctuations of the momentum. 

The relaxation of the time dependent flux (^) across the potential barrier to the time 
independent decay rate (|43|) is determined by the function 77 (t) in (44). In Fig. 2 the 
time dependence of r](t) is depicted for various temperatures. One sees that in the region 
of moderate damping the simple TST result Ttst = rr/(0) for the rate constant gives a 
satisfactory estimate of the true rate only for high temperatures. When the temperature 
is decreased rj(0) grows and depends strongly on the damping strength. Furthermore, for 
lower temperatures the average flux across the barrier becomes stationary faster for stronger 
damping. 



D. Relation to other rate formulas 

In the previous section we have calculated the probability flux across the potential barrier 
using the time dependent density matrix (|7]) with the initial preparation ([19|) . In particular, 
we have shown that the flux becomes time independent for times uu^t 3> 1 leading to the 
escape rate. Here, we want to regain the escape rate using rate formulas first introduced by 
Yamamoto [13] and Miller First, let us consider Yamamoto's rate formula 



r = lim -L rdX(Q[-q(-zX)}e[~q(t)]) p (68) 

t->oo tip JO 
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where the limit is understood as t ^> 1/cur. Here, the right hand side can be transformed 
to read 

-L dX(e[- q (-iX)}e[-q(t)})p = i([0[- g (t)],e[- g ]]) fl . (69) 



On the other hand, taking into account that 0(g) = 1 — Q(—q) one has from (48 



Im{C RL (t)} = -Im{C LL (t)} = ^ei-gWiehg])^. (70) 



Hence, we get from (BDj) 



r = ^feto{Qa.(«)}- (71) 



The result (|53|) can now be inserted into the above rate formula. First, from (|54]) and ( |35| ) 
one obtains for times u^f 1 

Im {log f | + 1=2 arctan [A(t)/S(t)] . (72) 



Thus, we obtain from (|53| 



*(«). 



tai m{CRl(t) }=^i^i_(n^ u „) (73) 



which combines with (|7T|) and the normalization ([41]) to yield the escape rate 

On the other hand, the time derivative Crl(£) given in ( [56]) determines Miller's rate 
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r= limC-RL^). (74) 



In the long time limit the imaginary part of Crl(^) becomes exponentially small and 

1 1 / 00 \ 

yields with ( [74]) again the rate (|43|). 

We note that for long times the flux-flux autocorrelation function ( |58"P becomes expo- 
nentially small. This indicates a constant flux across the barrier independent of the initial 
preparation of the nonequilibrium state in the metastable well. 
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V. CONCLUSIONS 



Within the path integral approach we have evaluated the time dependent density matrix 
of a metastable system in the vicinity of a barrier top when preparing the system at t = 



in thermal equilibrium on the left side of the barrier only (19). The explicit solution fl23|) 
is valid over a wide range of time excluding very long times and for high as well as for 
lower temperatures where quantum effects become important. The nonequilibrium state 
approaches an equilibrium state as one moves away from the barrier top. Condition fl3"B|) on 
the damping strength ensures that equilibrium is reached within the range of validity of the 
harmonic approximation for the barrier potential. 

In particular, we have studied the relaxation of the time dependent nonequilibrium state 
to the stationary flux state. We found that the corresponding time dependent normalized 
flux across the barrier is decaying in time. For very high temperatures the initial flux 
coincides with the transition state theory rate. For long times the flux coincides with the 
stationary decay rate of the metastable state which was shown to be identical with the well- 
known rate formula for thermally activated decay in the presence of quantum corrections. 
Furthermore, we have shown that the real time approach can also be used to evaluate 
correlation functions which are encountered in other rate formulas. 
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FIGURES 

FIG. 1. Width A(t) in position space of the nonequilibrium state for high temperatures 
hPtoo = 0.05 (thick lines) and lower temperatures h(3coo = 2.0 (thin lines) for a Drude model 
with wd/^o = 100. Solid lines indicate small damping with 7 = 0.1, dashed lines stronger damping 
with 7 = 3.0. 

FIG. 2. Time dependence of the average flux across the barrier rj(t) for two temperatures and 
a Drude model with ujd/loq = 100 and with various damping strengths. 
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